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Abstract. Perturbations due to round-off errors in computer modeling are discontinuous and 
therefore one cannot use results like KAM theory about smooth perturbations of twist maps. 
We elaborate a special approximation scheme to construct two smooth periodic on the angle 
perturbations of the twist map, bounding the discretized map from above and from below. 
Using the well known Moser’s theorem we prove the existence of invariant curves for these 
smooth approximations. As a result we are able to prove that any trajectory of the discretized 
twist map is eventually periodic. We discuss also some questions, concerning the application of 
the intersection property in Moser’s theorem and the generalization of our results for the twist 
map in Lobachevski plane. 

1 Introduction 

In recent years an increasing interest has been addressed to problems related to numerical simu¬ 
lations of complex dynamics. Classical mathematical results in this field provide the theoretical 
foundation for the numerical modeling of stable dynamical systems, while the analysis of math¬ 
ematical problems arising due to round-off errors was started relatively recently. Probably the 
first theoretical study was published in 1984 in Q. The results of this paper show that the 
influence of arbitrary small round-off errors can change the behavior of a complex (chaotic) 
system drastically, and that the resulting behavior depends very sensitively on the fine struc¬ 
ture of the phase space discretization corresponding to these errors. It is enough to mention 
here the period multiplication phenomenon and the localization under the action of such 
perturbations |j|]-|^], when trajectories that should normally be dense in the phase space remain 
confined to a small region. Starting from 1984 a large number of publications were devoted to 
the mathematical analysis of problems related to numerical simulations of complex dynamics 
(see, for example, @-0 and references therein). 

One of the most complex and intriguing problem arising in this field is the analysis of the 
influence of round-off errors for systems whose trajectories are neither stable nor unstable, or 
such that they have local components of this type. A typical example here is a rotation map 
F : (y3, r) —> {^ + a, r) of the two-dimensional plane IR^ around the origin. The study of 

the discretized version of this map was started in |^]. If the angle a G [0,1) is irrational then 
any trajectory of the map F is dense on a circle. However the action of round-offs disrupts 
these invariant circles, and a’priori it is not clear that there are no trajectories going to infinity. 
Moreover in Q an example of the rotation-like map with this property was constructed. Starting 
from that time at least 10 papers were dedicated to the discussion of this problem, however it 
was not solved. 
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Discretized plane rotation 


Figure 1: Invariant set (black dots) of the 0.0035-discretized pure plane rotation 


Numerical simulations (see Figure ||) show that the invariant set of the discretized rotation 
looks very much like invariant tori in the Kolmogorov - Arnold - Moser (KAM) theory. Therefore 
in 0 it was conjectured that in spite of the discontinuities of the perturbations due to the round¬ 
off errors some kind of KAM theory can be constructed for this problem, namely to prove that 
any trajectory of the discretized rotation is eventually periodic. In this paper we follow this idea 
and the main result is the following theorem. 

Definition 1.1 A map F of the annulus 

M(r_,r+) := {{(p,r) G IR^ : 0 < r- < r < < oo} 

into IR^ is said to be twist if in polar coordinates (0<<y9<l, 0<r< oo) it can be written as 

r ^ ^ + (mod 1) 

I _ y '■ '' 

where the function ‘h(r) is C^-smooth for a sufficiently large I {I = 6 is enough) and satisfies the 
inequality 

\d^{r)/dr\ > 0 . 

Notice, that in principle, the constant r_|_ might be equal to oo, in which case the annulus 
M(0, oo) coincides with the entire plane IR^. 

Theorem 1.1 For any constants 0 < r_ < r_ < r_|_ < f_|_ < oo there exists eo > 0 such that for 
each 0 < e < £0 if r^ < r < r^ then f_ < {D^ o F)'^{ip, r) < r+ for any n G 7L^ and any angles 
If G [0,1). Moreover, if we additionally assume that \d^{r)/dr\ > Const jr for any r > 0, then 
for each e > 0 any trajectory of the e-discretized twist map B is eventually periodic. 

Perturbations due to round-off errors in computer modeling are discontinuous and therefore 
one cannot use results like KAM theory about smooth perturbations of twist maps. We elaborate 
a special approximation scheme to construct two smooth periodic on ip perturbations of the twist 
map, bounding the discretized map from above and from below. Using the well known Moser’s 
theorem ||9| we prove the existence of invariant curves for these smooth approximations. The 
existence of the invariant curves yields by some special shadowing argument the boundedness 
of the trajectories of the discretized twist map. 


2 Twist maps with C'°°-smooth perturbations 

Consider a family of C^^-perturbed twist maps of the annulus ^(r_,r+) 

{ ip^ipF $(r) eh{ip, r, e) (mod 1) _ i) 

I V ■ 

depending on the small parameter e > 0. 

In this section we intensively use the following well known result due to Moser about the 
existence of invariant curves for the system (|2.lD. 





Figure 2: Counterexample for the intersection property. 

Theorem 2.1 ^ Let us assume that 

• \d^{r)/dr\ > 0 for any r- < r < r+, r+ — r_ > 1. 

• The functions ^,h,g G C°° with respect to all variables, and h and g are 1-periodie with 
respect to the angle variable (p. 

• For a sufficiently small eo for any |e| < eo any closed curve around the origin interseets 

with its image under the action of the map (intersection property). 

Then for any 0 < £ < Eq in the annulus A there exists an & -invariant elosed curve. 

Once again the assumption about the C'^-smoothness can be weaken to the C^-smoothness 
for a sufficiently large / (/ = 6 is enough). 

The last of the above assumptions is called the intersection property. Actually in the proof 
of this theorem a weaker condition is verified. Namely, it is enough to check up this condition 
only for closed curves {{<p,r) : r = r{Lp)} satisfying the inequality \dr/dip\ < (i.e., curves 

C^-close to circles). This follows from the estimate of the smoothness of the invariant curve 
and the fact that the intersection property is used in the proof only to check up that after the 
change of variables any circle intersects with its image (see p.2,15). Therefore we can use 
this weaker condition as a definition of the intersection property. 

Moreover it is straightforward to show that if any circle centered at the origin intersects with 
its image than any closed curve intersecting two concentric circles with the distance 2e between 
them also intersects with its image under the action of the above map. 

It seems that using this argument one can weaken the intersection property even more to the 
property that any circle centered at the origin intersects with its image. Moreover one might 
think that the intersection property for the circles yields this property for an arbitrary curve 
close to a circle. However, the following lemma shows that this is not true. 

Lemma 2.1 There exists a C°° map f of the annulus A, for whieh the interseetion property 
holds for any cirele centered at the origin, but there is a elosed curve arbitrary elose to a cirele 
which does not intersect with its image. 

Proof. Let the map / be a superposition of the rotation by 7r/2 around the origin and the 
linear map S, contracting with the coefficient 1 — a on the vertical coordinate (0 < a <C 1). 
Notice that we first apply the rotation. Clearly under the action of this map any circle centered 



at the origin intersects with its image at points with angles 0 and vr. Consider now a curve F, 
coinciding with a circle everywhere except small neighborhoods of points with angles 0 and vr. In 
these neighborhoods we consider very small outward perturbations of the circle (see Figure ^). 
The image of the curve F is the same as the image of the circle, except the neighborhoods of the 
points with the angles 7 r /2 and 37 r/ 2 . Clearly if our perturbation is less than o, then the image 
of the curve F does not intersect with the curve itself. ■ 


Notice that in the above example the intersections of the circles centered at the origin with 
their images occur at points with angles 0 and tt, and they are not transversal. If the intersection 
is transversal, then there exists a neighborhood of the circle such that for any closed curve in 
this neighborhood the intersection property holds. However this neighborhood is rather thin 


and it is simple to change the construction in the proof of Lemma 2.1 to make the intersections 
to be transversal by adding some expansion along the horizontal coordinate with the coefficient 
(1 + 6 ), provided a > 6 > 0. It is also straightforward to construct the twist map with the same 
property by the superposition of the map / with the rotation by the angle yr, providing the 
constant 7 is small enough. 

It is worth noting that the twist property is very essential for the existence of the invariant 
curve and one cannot extend above results for the case of the pure rotation, i.e., when <k(r) = 
Const. 


Lemma 2.2 Let <h(r) = a = Const. Then for any sufficiently small e > 0 there exist a pair of 
C°°-smooth functions h,g satisfying the above assumptions such that the perturbed system (2.1) 
does not have an invariant curve and moreover there are trajectories going to infinity. 


Proof. Consider the following perturbed family of maps, depending on the parameter e > 0 


ip ^ (f a -\- eh (mod 1 ) 
r —> r + eg{(p) 


( 2 . 2 ) 


where h G [0,1) is a constant and g{(p) is a C°° 1-periodic function. The system (2.2) is a skew 
product of two maps, and the first (angle) coordinate does not depend on the second (radial) 
one. Clearly we can choose the constant h in such a way that the number a + e/i will be 
rational. In this case the first map has a periodic trajectory (pi,ip 2 , ... ,(pk- Now setting g = 0 
we immediately see that all the trajectories are periodic (with the same period), and thus do 
not fill densely corresponding invariant circles. 

To show that at least some trajectories may go to infinity we need to define the function g 
in a bit more delicate way. We define the function g at the points (pi as +1 and extend it as a 
C^-function to the unit circle, such that it rapidly goes to —1 in small neighborhoods of the 
points ipi and is equal to —I otherwise. This definition clearly implies the intersection property, 
while the trajectory starting at the point ((/?i, 1 ) goes to infinity with the linear rate en. | 


3 Twist maps with discontinuous perturbations 

Definition 3.1 By an e-discretization, e > 0, of a set X C IR'^ we mean a choice of an ordered 
lattice (a collection of points) in the set X such that the distances between its neighboring 
elements do not exceed e. By an operator of e-discretization we mean a map : X — > X^ 
that associates to each point a: G X its nearest point on the lattice = D^{x) G X^ (if there 
are several such points we choose the point with the minimal index). Recall that the lattice 
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Figure 3: Approximation for the operator Di. (a) One-dimensional approximation, (b) Two- 
dimensional approximation. By thick (thin) lines we indicate the approximation from above 
(below). On figure (b) by thick black dots we indicate integer points, thick-lined small squares 
correspond to the images under the approximation from above of unit rectangles centered at 
integer points, without 25-narrow strips near their boundaries, indicated by thin lines. Thin 
small squares correspond to the images under the approximation from below. 


Ag is assumed to be ordered. The value of the parameter e > 0 is called the diameter of the 
e-discretization or the magnitude of the corresponding perturbation. 


Definition 3.2 By an e-discretized (perturbed) system for the map / with the phase space X 
we mean a pair (/g, Ag), where /g = Dg o /. 

Note that contrary to usual perturbation schemes, the discretized (perturbed) system here 
is given not on the original phase space A but on the discrete lattice Ag. 

In the case of round-off errors in computer modeling we deal with a uniform e-discretization, 
i.e., for e = 1/A, A G K+, the lattice Ag consists of all points x G A with rational coordinates 
with the denominator A. 

Denoting the twist map (B) by F we may represent the e-discretized twist map D^o F as 
follows; 

{ ip + $(r) -h l:h{p, r, e) (mod 1) 

I r —> r-|-e5((y?, r, e) ^ ’ 

The perturbations due to round-off errors that we consider are discontinuous, but there is an 
additional small parameter 1/r (for r S> 1) in the first relation, which arises due to the fact that 
the perturbation of the angle is inversely proportional to the distance from the origin. 

Our aim is to construct C°° smooth approximations for the discretization operator Dg. 
The hrst approximation is the approximation from above in the sense that for any point 
X G we have 

|D+(x)| > |Dg(x)|, 

while for the second approximation D~ - the approximation from below the opposite inequality 

\D~{x)\ < |Dg(x)| 

holds. Both of these approximations should satisfy also the property that their graphs lie in 
the 35e neighborhood of the graph of the operator Dg. Here 0 < 5 <C 1 is the parameter of 























the approximation. Notice that D^{x) = eDi{x/£) for any e > 0 and therefore it is enough to 
construct the approximations for e = 1. 

We define the the approximation from above in three steps. On the first step we define the 
following one-dimensional piecewise linear function, depending on a small parameter 0 < <5 1: 


ftit) ■■= 


2(5 4. 

l/2-2(5'' 


[-fH-t) 


if 0 < t < 1/2-2J 

i+ 2(5)+Z)i(t)+ 2(5 if -2(5<t-+>i(t)-i 

< —(5, t > 0 

5 + 5 ) + + 1 + (5 if — ^ — 6 < t — Di{t) 

< ^ — 2(5, t > ^ 
if f < 0 


The graph of this function is shown by thick lines on Figure ^(a). By the construction we have 


ftii) < Di{t) if t < 0 
iff >0. 


Thus |/i~(t)| > |Di(t)| for any t £ IR^. Straightforward estimates also show that the distance 
between the graphs of Di and /{*" is equal to 26. 

From the first sight it seems that the one-dimensional approximation from above for the 
discretization operator can be done in a much more simple way, by the function: 


r 2x if \x\ <1/2 

fi{t) := < X — 1/2 if X < —1/2 
[x + 1/2 ifx>l/2 


However for any twist map with the perturbation e/i(t/e) on each coordinate almost all trajec¬ 
tories go to infinity, and thus there is no invariant curves. 

On the second step we smooth out this piecewise linear function to obtain a new function 
having a given number of derivatives. The simplest way to do it is to apply the following integral 
smoothing operator: 

Scjh{t) \= f Sa{s)h{t — s) ds, 

J\s\<a- 

where 0 < cr << 1, and the kernel £ C°°, Sa{t) = 0 for |t| > a and 


t^Srrit) dt 


a for A; = 0 
0 for 0 < /c < Z ’ 


where I is a fixed number. It is shown (see details, for example, in |Q) that for any continuous 
function h the function S^jh £ C°° and its derivatives up to the order I are effectively estimated. 
Moreover we have a good control over the difference h — S^h, which enables to show that if the 
parameter a is sufficiently small, for example, cr = (5^, then applying this operator to //~ we 
obtain a function //" £ C°°, being the approximation of the discretization operator from above. 
The distance between the graphs will be somewhat enlarged, but is still less than 35. 

The multidimensional discretization operator Di acts on each coordinate independently of 
others. Thus on the last step of the approximation procedure we apply the function //" on 
each coordinate independently, which defines the smooth define approximation Df for the two- 
dimensional discretization operator Di. The action of this two-dimensional transformation, 
shown on Figure |3|(b), is the following. First, it maps almost linearly unit squares centered 
at integer points (except 25-narrow strips near their boundaries, indicated by thin lines on the 
figure) to small squares of size 5, located near the corresponding integer points and shifted 





outward with respect to the origin. Second, the mentioned above strips are expanded to the 
complement of these small squares. Since /{*" is the approximation from above for the one¬ 
dimensional discretization operator, similar property holds also for : 

|Z)+(x)| > \D^{x)\ 

for any x G H^. Thus, according to our definition, is the approximation from above for Di. 
Straightforward estimates also show that the distance between the graphs Di and Df is not 
more than to 3(5. 

In a similar way one can construct the approximation from below ■ The corresponding 
one-dimensional version is defined as 


^ t 

1/2+5'' 



if 0 < t < 1/2(5 




if 5 < t i < 2(5 

tG(* 

-Dlit) - i + 2(5 

)+Dlit)-25 

if — ^ + 26 < t — Dlit) 




< i + 5, t i 

¥(* 

1 

1 

1 

if 5 < t — Dlit) + ^ < 2(5, 


t > \ 
ift<0 

The graph of this function is shown by thin lines on Figure ^(a). By the construction we have 

flit) > Di{t) if t < 0 
flit) < Diif) if t > 0 

for any point \t\ > 1. Thus \fiit)\ < |Tli(t)| for |t| > 1, which makes f^ the approximation 
from below everywhere except the central interval [—1,1]. Notice that there is no continuous 
one-dimensional approximation from below with nondegenerate derivatives for Di for all t G . 

The second and the third steps of the procedure are exactly the same as in the construction 
of the approximation from above. The scheme of the action of the approximation from below is 
shown by thin-lined small squares on Figure ^(b). It is worth noting that the operator does 
not satisfy our requirement to be the approximation from below in the 1 -|- 2(5-square centered 
at the origin. However we shall apply our approximation scheme only for the annulus A with 
nonzero smaller radius. Therefore we do not need to define the discretization operator near the 
origin. Another problem is to verify the approximation from below property for unit squares 
intersecting with the coordinate axes. The operator restricted to such squares enlarges one 
of the coordinates and decreases another one. However, take an integer point (0, n) and notice 
that 

(n — (5)^ + 6"^ = — 2(5(1 — 6) < n^. 

Therefore the emphasized small square lies entirely inside of the circle of radius n centered at 
the origin. Thus the effect of the decreasing is stronger, which proves the desired property. 

The proof of the intersection property for the approximation from above is as follows. Con¬ 
sider a circle centered at the origin with a sufficiently large radius. Then locally this circle 
is arbitrary close the the straight line. Take the segment of the circle with the “slope” — 1. 
Clearly there is a point i^,r) in the lattice of “half-integer points” -|- (1/2,1/2)) arbitrary 
close to this segment in the same “half-plane” as the origin. Then by the construction of our 
approximation there is an intersection of the image of the circle with the circle itself near the 
image of the point i<f,r). Indeed in the 2e(5-neighborhood of the point i(p,r) the function Df 





takes both values {ip-,r — e/2) and {(p+,r + e/2), which yields the intersection. In the same 
way we can deal with a curve sufficiently C^-close to the circle. In the same way one can prove 
the intersection property for the approximation from below D~. 

Therefore the both maps Df o F satisfy the conditions of Theorem |2.l| , and thus they have 
invariant curves, defined by the relations r = rj-((^). In spite of the fact that our approximation 
is the approximation from above (from below), the angles corresponding to the points DfoF{x) 
and o F{x) are different. Therefore we cannot apply the idea that every time the trajectory 
of the approximating map is further from the origin or at least not closer (on closer in the case 
of the approximation from below) than the discretized one, being at the same time uniformly 
bounded (by the invariant curves). 

To prove the boundedness we need an additional trick. For any point x G TZ? consider its 
“shadow” point x^ on the invariant curve (the point with the same angle). Now comparing the 
images of these points o F{x^) and o F{x) we obtain the boundedness of the trajectory 
of the discretized map. 

This completes the proof of the first part of Theorem 0. 


To prove the second part we restrict our twist map F to the annulus ^(r_,r+) such that 
1 <C < r_ < 2R < r_|_ < with a sufficiently large R. Then, after the following change of 
variables: {'-p,r) —> ((/?, r_r), we obtain the twist map of the new annulus .4(1, r_|_/r_) with the 
perturbation of order e/r_, which can be made arbitrary small by the choise of the constant r_. 
Therefore the argument in the first part of the proof yields the boundedness of any trajectory 
of the discretized map with arbitrary e < oo. On the other hand, in a bounded region there is 
only a finite number of points of the discretized lattice, which yields the eventual periodicity. 
Theorem 1.1 is proven. | 


The similar construction can be applied for the “integer part” discretization, when each 
point is mapped to its integer part. In spite of of some differences, this type of discretizations 
leads to the same eventually periodic behavior. On the other hand, the truncation operator for 
arbitrary small truncation error leads to quite different behavior: almost any trajectory of the 
truncated system goes to zero with time |^]. 


4 Discretization of the twist map in Lobachevski’s metric 

In this section we discuss some our findings about the properties of discretizations of the twist 
map in Lobachevski’s metric. The Lobachevski plane IL 2 can be realized as the upper semi-plane 
y > 0 of the Eucledian two-dimensional plane IR^ with the metric = {{dx)‘^ + {dy)'^)/y^, which 
we denote by p(-, •). We shall consider two different discretizations of IL 2 : one with respect to the 
usual integer lattice and another one with respect to TL^ = {{x,y),x G y = 2^, A: G ZL}. 
The reason for the introduction of the second lattice is that it is uniform with respect to the 
Lobachevski’s metric, which is not the case for the lattice Notice that the lattice ZL^is not 
locally finite. 

Observe that the distance in the Lobachevski’s metric between any two horizontal straight 
lines vanishes as x ^ ± 00 . Therefore any vertical straight line plays a role of the radius 
from the infinity y = 00 , while any horizontal straight line is a circle centered at the infinity. 
Consider a map F : IL 2 —> IL 2 defined as F{x,y) = (x -|- <h(y),y). The discretization operator 
F : IL 2 —> U we define as the closest point map in the Lobachevski’s metric. 

Proposition 4.1 The intersection of ZI? with any horizontal line with the integer coordinate 
y ^ ZZ is invariant with respect to the iZL^ U ZLjf)-discretized map DoF. 



Proof. It suffices to show that D{x, n) = (m, re) for any numbers x G H and n & and some 
integer m € 7L. Let x = m + S, where m € Th and |(5| < 1/2. Then for D{x,n) = {x',y') we 
have: 

p{{x,n), {x',y')) > 

Therefore if y' ^ n and re > 1, then 

p{{x,n),{x,y')) > |ln((re +l)/re)| > - - > —• 

re 2re^ re 

On the other hand, p{{x,n), {x',y')) < \S\/n for {x',y') = {m,n). 

It remains to check the case re = 1 and y' = 2~^ with k > 1, k G 7Zi. In this case 


fy' ^ 
In y 


= \ ln{y'/n)\. 


p{{x,l),{x',y')) > 



ln2 > \6\ > p((x, l),(rre,l)). 


It is of interest that the intersection of the second lattice with horizontal lines is not 
invariant with respect to the discretized map. The proof of this statement and the study of the 
discretization of finite radius rotation maps in Lobachevski’s metric is technically sufficiently 
complex and will be published elsewhere. 
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